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Abstract 

We present an estimate of certain higher order corrections to the contribution of 
the charm triangle loop in the inclusive B — >> X s ^ decay rate, recently discussed by 
Voloshin. We find that these corrections are minute and hence the result found by 
Voloshin, although small, is quite robust. 



The concepts of heavy quark universality, symmetry and effective field theory have been 
exhaustively applied to inclusive B decays [lj . By these means, the inclusive B —>■ X c lv and 
B — > X s 7 decays can be related to the underlying (perturbatively computable) b — > civ and 
b —* sj quark decays with only 0(Aq CD /77t|) soft physics and 6-quark binding corrections. 
These results have been brought into question by a recent paper of Voloshin [[|, where a 
non-perturbative correction to the inclusive rate for B — > Xgj is identified which scales as 

AqcdMc- 

The appearance of non-perturbative corrections apparently missed in the heavy quark 
effective theory (HQET) treatment of inclusive B decays is unsettling and deserves a more 
thorough discussion. In this note we show how, in the unphysical limit where m 2 3> ^/i^Aqcd, 
it is possible to understand how the Voloshin correction arises in the context of HQET. In 
this limit, in fact, this correction is formally less than terms of 0(m 2 /m 2 ) and hence is 
under control. In real life, however, m 2 ~ m;AQCD and, in principle, Voloshin's correction 
is subject to considerable uncertainty. In practice, however, we show that all corrections to 
Voloshin's result are quite negligible, so that indeed his computation provides the dominant 
contribution. 

The corrections of 0(AQ CD /ml) identified by Voloshin [fj] arose by considering the con- 
tribution of the gluon-photon penguin graph, shown in Fig. 1. After applying a Fierz 
transformation to the SiTliiF-LCLri^bL four-quark operator in the underlying effective Hamil- 
tonian this graph is proportional to the famous AVV triangle diagram. Because of the GIM 
mechanism, there is no anomaly in the axial vertex, and the result is necessarily non-singular 
in (fc 7 + kg) 2 . In the limit of vanishing gluon momentum considered by Voloshin 0, this 
contribution must scale as V^V^/m 2 for the up-type quarks i running in the loop. In his cal- 
culation, Voloshin ignores the up quark loop because of its miniscule CKM matrix elements 
and because, in no sense, can m u be considered larger than the typical gluon momentum in 
the problem. Since V* s V c b ~ —V t * s V t b and m 2 3> m 2 c) the c quark loop in Fig. 1 dominates 
and, indeed, the dominant effective operator for the b — > sjg process scales as 1/m 2 . This 
is the origin of Voloshin's result for the effective Lagrangian for this process: 

47i S79 = ^V2G F V^V cb (s L ^b L )^G?d x P,„. (1) 

Other observations regarding effective Lagrangians of this type have been made previously. 
The effect on the exclusive B — > K*-f decay has been studied in Ref. ||, while the authors 
of Ref. have studied the effect of the b — > s^g process on the photon energy spectrum. 

The interference of the above term with the leading perturbative contribution for the 
b — ► s'-f process, given by the effective Lagrangian || 

A_ 7 = 6 ^ V^CMmbS^MF^, (2) 

lbir z a/2 

in the operator product expansion (OPE) for the inclusive rate yields an effective amplitude 
for the b — > sj process 



71-7 = -^G 2 F ml±Re[V; s V cb V ts V t lC 7 ] 
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Figure 1: Triangle diagram for the b — > s^g process. We have omitted a second diagram 
with the gluon and photon interchanged. 



This produces a correction to the B —>■ Xgj rate: 

5T(B^X sl ) 1 



^ 



r(S^X s7 ) 27C 7 m 2 c 



-0.025. 



(4) 



The numerical result above follows by using for the coefficient, Cj, of the leading operator 
mediating the B — > X s j transition, C-j ~ —0.3 ||, and using the standard evaluation for the 
strength of the chromomagnetic interaction of the 6-quark inside the B hadron 0, 



A*! = \{B\bga rv Gl^b\B) = 2(M 
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0.4 GeV 2 . 
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Although (H) is a small correction to the inclusive rate, its sensitivity to the scale m 2 c 
brings into question the expectations from HQET. To understand what is going on, it is useful 
to examine more fully the gluon-photon penguin graph considered by Voloshin. Because 
the AVV graph has been analyzed in detail by Adler ||, it is straightforward to consider 
corrections to the effective Lagrangian ([[]) or, equivalently, the amplitude (H). 

The triangle graph in question is proportional to the tensor 



lfiaf3 — A^p kg) {^^ii/afiykg ' kj)(k^ Kg) [^i/fMarkj p 

where the invariant function, I(k y ,k g ), is given by 
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with rrii being the mass of the quark in the loop. In the limit where the gluon 4-momentum 
vanishes, Eq. ([?]) reduces to the 1/m 2 factor alluded to earlier. In fact, the gluon 4-momentum 
is never vanishing. It is of order of the typical momentum of the B meson constituents, which 



is of C(Aqcd)-Q Obviously, for the w-quark loop, it makes no sense to consider the k£ — > 
limit and so the scale l/m 2 ^ never enters the problem. Given the very small magnitude of 
Vus^ub relative to V* s V c b, it is perfectly sensible — as Voloshin does — to neglect the w-quark 
loop altogether. 

The situation is also clear for the t-quark loop. In this case m t is really much larger than 
all other scales. Given that the k^ — > limit is appropriate, the photon-gluon penguin graph 
involving the t-loop will contribute an effective interaction like that of Eq. ([I]), but scaled 
by l/m 2 . This interaction, however, gives a negligibly small correction to the B — > X s 7 
rate. For the c-quark loop, however, the effective k*t — > limit which yields the b — > s'yg 
Lagrangian of Eq. (JI]) does not appear so safe. 

It is certainly possible to drop the k 2 x(l — x) term in the denominator of Eq. (^) relative 
to ml, since m 2 3> Aq CD . However, in the B rest frame |k 7 | ~ rrib/2, thus the 2xyk g ■ fc 7 
term in Eq. (|7|) is of C^m^AqcD), which is not really small compared to m 2 . Voloshin's result 
(Eq. ([!])) neglects this term altogether. 

Before evaluating corrections to Voloshin's formula, it is useful to make some general 
remarks. The approximation of also dropping the 2xyk g ■ kj terms in Eq. ([?]) for the c-quark 
loop is only tenable in a world where m 2 ^ Aqcd^^V If this were the case, the Voloshin 
correction to the B — > X s j inclusive rate, given in Eq. (f|), would not violate, per se, the 
heavy quark expansion, since 

^ - *f7) „ i£ „ ^ocr, <<; K (8) 

T(B^X sl ) m 2 m 2 m 2 ' K ' 

On the other hand, dropping the 2xyk g ■ /c 7 term is never a good approximation in the 
heavy m^ limit (m^ — > oo, with m c fixed) envisaged in the HQET. In this limit, the effective 
Lagrangian ([I]) scaled by l/m 2 never arises, since the proper limit for Eq. (|?p is not l/m 2 but 
—6/kg ■ fc 7 . Indeed, in this limit the photon-gluon penguin contribution vanishes identically 
because of the GIM mechanism, V*Vib = 0. 

A class of systematic corrections to Voloshin's result can be computed by retaining the full 
function 7(/c 7 , k g ) in the photon-gluon penguin graph. To be more precise, since m 2 c 3> Aq CD , 
it suffices to consider 

rl rl — x j 

J(fc 7 , k g ) = 24 / xdx / ydy— — - — — . (9) 

Jo Jo [m 2 c — 2xyk^ ■ k g \ 

Expanding the denominator in powers of fc 7 • k g /m 2 identifies a progressive set of higher di- 
mensional operators which contribute to the h — > s^g Lagrangian. This expansion effectively 
replaces the G v a x d\F^ v /m 2 operator in Eq. (|1|) by 

1 ~ 1 

— ~G V id\Fn U — > — -G u id\Fn V 
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+ ^(id a idpid x F, v ){iD a iD l3 G: x ) + ... (10) 



1 To be more precise, as we will argue later, |k g | ~ O(Aqcd)- 





Figure 2: The cut diagram which, together with its complex conjugate, yields the contri- 



bution of C 



b-ts-yg 



to the inclusive B — > X s ^ rate. 



The interference of these higher dimensional terms with the leading order operator (|j) for 
the b — > S7 process in the operator product expansion formula for the inclusive rate provides 
the desired set of corrections to the Voloshin result. 

In practice, it is simpler to compute these corrections by first calculating the correlator of 
the leading order contribution to b — > sj, Eq. (Q), with the full c-quark photon-gluon penguin 
graph and then expanding the result in powers of the gluon 4-momentum. A tedious but 
straightforward calculation of the correlator depicted in Fig. 2 yields the following formula 
for the effective amplitude for the b — > sj process: 
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Here J a \ involves an integral of I(k^, k g ) over the photon and s-quark phase space and is 
given by 
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where P is the b quark 4-momentum. 

The tensor J a \ can be expanded in terms of scalar functions which depend on the invari- 
ants which are left over after the phase space integrations (kg, P ■ k g , and P 2 ): 



JaX = JlTjaX + ^2-Pq-Pa + Jzk ga k g \ + J±{k ga P\ + k gX P a ). 



(13) 



The scalar functions J« are easily identified by contracting J a \ with the gluon and/or the b 
quark 4-momentum. They involve combinations of the phase space integrals of I(/c 7 , k g ) , 
(k-y ■ fc s )J(fc 7 , kg) , and (fc 7 • /c 9 ) 2 /(A; 7 , k g ). These integrals are readily computed in a power 
series in the gluon momentum. We will illustrate this with the phase space integral of 
J(fc 7 , k g ), which we perform in the b rest frame: 
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The gluon energy in the above is typically much smaller than the gluon momentum. Roughly 
speaking, the B matrix elements involving the gluon energy are of order the difference in 
kinetic energy of the b quark before and after absorbing a soft gluon. Hence the gluon 
energy is of order E g ~ AT b ~ AQ CD /m b , while |k ff | ~ AqcD~the gluons are predominantly 
spacelike. Dropping the gluon energy E g relative to |k s | and expanding the logarithm in 
powers of m b AqcD/ml yields a rapidly converging power series for (I(/e 7 , k g ))\ 
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The higher order terms in the series involve integrals of x n y n , which fall off asymptotically 
as v/?/n 3 / 2 4 n+1 . 

Using the decomposition (j^) in Eq. ( |TT] ) and retaining only the terms that do not vanish 
by the equations of motion gives for the b — > s'y amplitude the expansion 
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Calculations analogous to those that led to Eq. ([TJ) lead to the following expressions for the 
leading behavior of the scalar functions J\, J%, and J4 in the b rest frame Q: 
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Although the J4 term in Eq. (|16D is nominally linear in the gluon 4-momentum, it actually 
gives a correction of 0(AQ CD /ml) upon using Faraday's law for the gluon fields (k x E a = 
£ 5 B a ), since E g ~ AQ CD /m6. Thus the leading correction to Voloshin's result for T b ^, sl 
given in Eq. (3) arises only from the J\ and J3 terms and involves operators quadratic in 
the gluon 4-momentum. Using Eq. fli~7|) , one identifies this correction as 
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2 The J2 term drops out because of the equations of motion. 



Although one cannot relate the above matrix element to a property of the B mesons, as was 
done for the leading matrix element, we expect it to be of 0(Aq CD ). However, because of 
the tiny numerical coefficient, even though ?7?|Aq CD ~ m^, the above correction is totally 
negligible. So Voloshin's result @ is robust, at least as far as these corrections go. 

There are, of course, many other corrections to Voloshin's result. However, these are all 
of 0(AQ CD /m*), down by a factor of Is^ CTl /m 2 c relative to the leading term calculated by 
Voloshin. Some of these arise by considering the full photon-gluon penguin graph, e.g. the 
J4 term in Eq. fll8|). Others come from terms where £ V ~ > 1 ' S79 is correlated with itself in the 
OPE for the total width. Yet others come from photon-gluon penguin graphs in which two 
or more soft gluons are emitted. Although all of these contributions involve unknown matrix 
elements, because Aq CD <C m 2 c it is reasonable to expect that they also should yield quite 
small corrections to Eq. (f|). 

Our analysis gives some confidence that the 1/m^ nonperturbative corrections to the 
B —* Xs'-f process calculated by Voloshin || are an accurate estimate of these corrections. 
Unfortunately, these interesting corrections are numerically small for this particular decay. 
Their existence, however, raises the interesting question of whether other 1/rajj? corrections 
may give more sizable contributions to other processes. We hope to return to this issue at a 
later date. 

Acknowledgments 

We are grateful to Mark Wise for useful discussions on these matters. He and his col- 
laborators have arrived at similar conclusions to ours. [See Z. Ligeti, L. Randall, and M. 
Wise, CALT-68-2097, MIT-CTP-2611, |hep-ph/9702322| .1 We thank them for discussing their 
results prior to publication. We also would like to thank Misha Voloshin for a helpful dis- 
cussion. One of us (S.N.) is grateful for the hospitality of the UCLA department of Physics 
and Astronomy, where part of this work was completed. This work was supported in part 
by Department of Energy grant DOE-FG03-91ER40662, Task C. 

References 

[1] Nonperturbative contributions to the B — > X s 7 rate have been discussed in the context 
of HQET by A. F. Falk, M. Luke, and M. J. Savage, Phys. Rev. D49 3367 (1994). For 
reviews of heavy quark effective theory, see, e.g. M. B. Wise, in Gatlinburg 1994, Pro- 
ceedings, Radiative Corrections ( |hep-ph/9411264D , T. Mannel, talk given at the Work- 
shop on Heavy Quarks at Fixed Target (HQ 96), St. Goar, Germany, Oct. 3-6 1996 
(|hep-ph/9611411|) . 



[2] M. B. Voloshin, University of Minnesota preprint TPI-MINN-96/30-T flhep- 
ph/9612483D . 



[3] A. Khodjamirian, R. Ruckl, G. Stoll, and D. Wyler, preprint WUE-ITP-97-001, MPI- 
PhT/97-5, Zurich Uni-31/96, Technion-ph-96-24, |hep-ph/97023T^ . 

[4] A. Ali and C. Greub, Phys. Lett. B361 46 (1995). 



[5] B. Grinstein, R. Springer and M. Wise, Nucl. Phys. B339 269 (1990). 

[6] K. Chetyrkin, M. Misiak and M. Miinz, preprint ZU-TH 24/96, TUM-HEP-263/96, 
MPI/PhT/96-123 ( |hep-ph/96123T3| ). 

[7] See, e.g. A. F. Falk and M. Neubert, Phys. Rev. D47 2965 (1993). 

[8] S. L. Adler, Phys. Rev. 177 2426 (1969). 



